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ABSTRACT 

We present Gran Telescopio Canarias (GTC) optical transit narrow-band photometry of the hot-Jupiter exoplanet 
XO-2b using the OSIRIS instrument. This unique instrument has the capabilities to deliver high cadence narrow-band 
photometric lightcurves, allowing us to probe the atmospheric composition of hot Jupiters from the ground. The observa- 
tions were taken during three transit events which cover four wavelengths at spectral resolutions near 500, necessary for 
observing atmospheric features, and have near-photon limited sub-mmag precisions. Precision narrow-band photometry 
on a large aperture telescope allows for atmospheric transmission spectral features to be observed for exoplanets around 
much fainter stars than those of the well studied targets HD 209458b and HD 189733b, providing access to the majority 
of known transiting planets. For XO-2b, we measure planet-to-star radius contrasts of R P !/R*=0.10508±0. 00052 at 
6792 A, 0.10640±0.00058 at 7582 A, and 0.10686±0.00060 at 7664.9 A, and 0.10362±0.00051 at 8839 A. These mea- 
surements reveal significant spectral features at two wavelengths, with an absorption level of 0.067±0.016% at 7664.9 
A due to atmospheric potassium in the line core (a 4.1-<r significance level), and an absorption level of 0.058±0.016% 
at 7582 A, (a 3.6-ct significance level). When comparing our measurements to hot- Jupiter atmospheric models, we find 
good agreement with models which are dominated in the optical by alkali metals. This is the first evidence for potassium 
in an extrasolar planet, an element that has long been theorized along with sodium to be a dominant source of opacity 
at optical wavelengths for hot Jupiters. 
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1. Introduction 

Transiting hot-Jupiter exoplanets provide an excellent op- 
portunity to detect and characterize exoplanetary atmo- 
spheres. During a transit event, both the opaque body of 
the planet and its atmosphere block light from the parent 
star, allowing the radius and its wavelength dependence 
to be accurately measured, giving rise to the detection of 
atomic and molecular features in the atmosphere. In addi- 



* Based on observations made with the Gran Telescopio 
Canarias (GTC), installed in the Spanish Observatorio del 
Roque de los Muchachos of the Instituto de Astroffsica de 
Canarias, in the island of La Palma, and part of the large ESO 
program 182.C-2018. 



tion, secondary eclipse measurements can be used to obtain 
an exoplanet's emission spectra, where such properties as 
the temperature, thermal structure, and composition can 



be measured ( 


e.g.|Deming et al. 2005 


Grillmair et al. 2007 


Charbonneau et al. 2008 ; Rogers et a 


2009;|aing& Lopez- 


Morales 


2009 


) while orbital phase curves can probe the 



For highly irradiated planets, the atmosphere at opti- 
cal wavelengths represents the "window" into the planet 
where the bulk of the stellar flux is deposited. Thus, the 
atmospheric properties at these wavelengths are directly 
linked to important global properties such as atmospheric 
circulation, thermal inversion layers, and inflated planetary 
radii. All of these features may be linked within the hot- 
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Jupiter exoplanetary class ( Fortney et al.|2008 ), and trans- 
mission spectra can provide direct observational constraints 
and fundamental measurements of the atmospheric proper- 
ties at these crucial wavelengths. In the optical, hot- Jupiter 
transmission spectra are thought to be largely dominated 
by the opacity sources of the alkali metals sodium and 
potassium, with both species predicted to be present since 
the initial model calculations of highly irradiated gas giant 
exoplanets ( Seager fc Sasselov|2000 Brown|200T Hubbard 



et al.pOOl ) . Other strong but still unidentified optical ab 



sorbers could also play a role in the formation of strato- 



spheres ( |Burrows et al.|2007 



^ |Knutson et al7[2008| ) , with the 
leading candidates being "TiO/ VO and sulfur compounds 



(IHubeny et al |2003[ [Fortney et al.|2008||Desert et al. 12008, 

Z ahnle et al. 2009| ). No candidate is yet a clear favorite 
to universally produce hot- Jupiter stratospheres, both the 



TiO and su lfur explanations have unresolved issues ( Spiegel 
et al.||2009[) and stellar act ivity levels could also be impor- 
tant |Knutson et al.|p20To] ) 



To date, only two exoplanets have been ideal for op- 
tical atmospheric studies, HD 209458b and HD 189733b, 
mainly due to the stellar brightness (V~ 7.8) along with the 
large transit signals and low surface gravities. Atmospheric 
sodium has been detected in HD209458b from data with 



HST STIS ( |Charbonneau et alj[2002j |Sing et al.||2008a M 
and confirm ed by ground-based spectrographs ( |Snellen| 
et al.||2008[ ). In HD 189733b, a ground-based detect ion of 
sodium has also been made by Redfield et al. (2008), with 
the wider ove rall optical spectrum consistent with a high 



altitude haze (Pont et al. 2008 Sing et al. 2009 Desert et al 



2009), thought to be due to Rayleigh scattering by smal l 



condensate particles (Lecavelier des Etangs et al. [2008) 
While sodium has been detected in both of these planets, 
the other important alkali metal potassium has not, though 
in principle both the ground-based and HST spectra should 
have sufficient sensitivity. 

Here we present the first results from GTC OSIRIS for 
the exoplanet XO-2b, part of our larger spectrophotomet- 
ric optical survey of transiting hot Jupiters. XO-2b is a 
0.996i?j up , 0.565 Mj up planet in a 2.6 day orbital period 
around a n on-active V=11.2 early K dwarf star ( |Burke 



et al. 2007). XO-2b falls in the transition zone between 
the proposed "pM/pL " hot- Jupiter classification system of 
Fortney et al. (20081, and has been recently r eported to 



et al. 


2009) 


mospl 


leric i 



Our overall program goals are to detect at 



mosphcrcs, enabling comparisons between the features ob- 
served (comparative exoplanetology) . OSIRIS is an optical 
imager and spectrograph on the new 10.4 meter GranTeCan 
(GTC) telescope, capable of extremely high-precision differ- 
ential narrowb and fast photomet ry, even for fainter transit 
hosting stars ( |Col6n et al.||2010| . In our program, transit 
events are observed with a narrowband tunable filter, ca- 
pable of simultaneously providing both the spectral resolu- 
tion and extremely high photometric precision necessary to 
detect atmospheric transmission spectral features. Multiple 
wavelengths are chosen to be maximally sensitive to the ex- 
pected hot- Jupiter atmospheric features, including sodium, 
potassium, and titanium oxide. An important advantage 
to spectrophotometry is that the absolute transit depths 
of each wavelength are measured, allowing for atmospheric 
features to be compared together alongside existing broad- 
band data at other wavelengths. In addition, transit timing 




Fig. 1. GTC Osiris CCD frame of XO-2 from 4 Dec 2009 
taken with the tunable filter, with both CCD1 (left) and 
CCD2 (right) shown. XO-2 is the circled bright star lo- 
cated near the center of CCD1, with XO-2's stellar com- 
panion, used for our reference star, nearby and directly 
overhead. The z-scale of the image has been set to best 
enhance the weak emission sky lines, which can be seen as 
"rings" around the image. Vignetting from components in- 
side the instrument produce the artifacts seen at the very 
left edge of the image. 



variation studies can also be made along with refining an 
exoplanet's system parameters. In this paper, we describe 
the GTC XO-2b observations in §2, present the analysis of 
the transit light curves in §3, discuss the results and com- 
pare to model atmospheres in §4, and conclude in §5. 



2. Observations 

We observed with the GTC 10.4 meter telescope installed 
in the Spanish Observatorio del Roque de los Muchachos 
of the Instituto de Astrofi'sica de Canarias on the island of 
La Palma, using the Optical System for Imaging and low 
Resolution Integrated Spectroscopy (OSIRIS) instrument 
during three separate nights in service mode. 

2.1. GTC OSIRIS instrument setup 

The GTC OSIRIS instrument ( |Cepa et al.|2000| [2003] ) con- 
sists of a mosaic of two Marconi CCD detectors, each with 
2048x4096 pixels and a total FOV of 7.8'x8.5', giving a 
plate scale of 0.127"/pixel. The CCD has a pixel well depth 
of ^100,000 electrons and a 16 bit A/D converter, which 
saturates at 65,536 counts. For our observations, we choose 
the lxl binning mode with a full frame readout at a speed 
of 500kHz, which produces a readout overhead time of 31.15 
seconds. The readout speed of 500kHz has a gain of 1.46 
e~ / ADU and a readout noise of 8 e~ . While 500kHz is the 
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Table 1. GTC OSIRIS XO-2 observations 



HJD start 


A (A) 


Number of 


Integration time 






exposures 


(sec) 


2455149.546 


6792.0 


249 


20.47 


2455170.452 


7664.9 


147 


25.47 


2455170.454 


7582.0 


149 


25.47 


2455259.411 


8839.0 


223 


40.67 



noisiest readout speed (100 and 200kHz are available), it 
allows for the most counts per image to be obtained due to 
a higher gain. Our initial service mode observations have 
been taken with the initially available full-frame readout 
mode, a conservative observing approach for our initial ob- 
servations, with highly efficient windowing now becoming 
available. In future observations, windowing (sub-arrays) 
will eliminate nearly all CCD readout time and provide 
better time-sampling with a near 100% duty cycle. 

We used the tunable filter (TF), which allows for 
narrow-band imaging capable of specifying the wave- 
length^) of our target and reference stars. The TF consists 
of a Fabry-Perot etalon, whose cavity separation can be ad- 
justed, giving rise to a wide range of resolving powers, with 
OSIRIS capable of R=300 to 1000. The unwanted orders 
from the etalon are then blocked by a set of broadband and 
mediumband (100-600 A FWHM) filters. Different path- 
lengths through the optical system produce a wavelength 
dependance within the field of view, such that each im- 
age has "rings" of constant wavelength with respect to the 
TF-center, located near the middle of the detector system 
(see Fig.JT]). To provide a maximum spectral resolution and 
largest contrast for atmospheric signatures, we choose to 
observe with the narrowest width available to all wave- 
lengths of 12 A. 

For the XO-2b observations, we placed both the target 
and a chosen nearby reference star, XO-2's stellar compan- 
ion XO-2B, at the same radial distance from the TF-center 
such that they would be observed at the same wavelength 
(see Fig. [I]). The companions XO-2 and XO-2B are sepa- 
rated by 31" and cataloged as each having the same spectral 
type (K0V) and have very similar magnitudes and colors, 
making XO-2B an ideal reference star. When searching for 
the strong atmospheric potassium line, successive images 
were taken alternating between two different wavelength 
positions, carefully chosen to also avoid bright sky lines and 
strong telluric absorption lines, yet located close enough 
together in wavelength to be within the same order-sorter 
filter. As the tunable filter can scan between wavelengths 
within a given order-sorter filter very rapidly (~1 ms) and 
accurately (~0.1 A), we were able to observe at two potas- 
sium line wavelengths during a single transit while avoid- 
ing the typical long overheads associated with switching 
between traditional filters. 



2.2. Observing Log 

XO-2 was observed during the three transit events of 13 
November 2009, 04 December 2009, and 03 March 2010 
(see Table [TJ. 

13 Nov. 2009: The seeing was very good during the ob- 
servations, with values ranging from 0.6" to 0.9". The filter 
was tuned to a sing le wavelength of 6792.0 A, with 20.467 
second exposure times used during the observing sequence 



with about 250 exposures acquired. During the observa- 
tions, there was an interruption of 20 min due to problems 
with the M2 collimator, causing ingress to be largely lost. 
The out of transit baseline flux was also truncated, com- 
pared to the other nights, as time was lost with both the 
guiding setup and the TF entering a forbidden range at the 
end of the night, where calibration becomes very unreliable. 
However, we found there was still ample time to accurately 
measure the baseline flux both before and after the transit, 
and the interruptions and delays did not limit our analysis. 
During the observations, the airmass decreased during the 
night from values of 1.678 to 1.08 

4 Dec 2009: For the transit, XO-2 was observed at 
two TF wavelengths, switching successive images between 
7664.89 and 7582 A. The night was photometric, weather 
conditions very stable, and the seeing was very good dur- 
ing the sequence, with values ranging from 0.7" to 0.9". 
Approximately 150 exposures at each wavelength were ac- 
quired. The total exposure time was 25.466 seconds for 
each image, producing successive frames at the same wave- 
length every 113.2 seconds, and an overall duty cycle of 
45%. During the observations, the airmass decreased dur- 
ing the night from values of 2.158 to 1.081. 

3 March 2010: Weather conditions were very stable, and 
the seeing was good during the sequence, with values rang- 
ing from 0.7" to 0.9". The filter was tuned to a single wave- 
length of 8839. OA, with 40.67 second exposure times used 
during the observing sequence with about 230 exposures 
acquired. During the observations, the airmass increased 
during the night from values of 1.078 to 1.639. 

For all three nights, exposure times were not modified 
during the sequence, maintaining the peak counts around 
30000-35000 ADUs for both the reference and target stars. 
When needed, defocusing was used to avoid saturation in 
improved seeing conditions. No recalibration of the TF was 
needed during the observing sequences, as the field rotator 
moved in a range where no variations are expected. From 
data obtained during TF commissioning, possible variations 
in the central wavelength of the TF in the present observa- 
tions are smaller than lA, with real variations likely even 
lower. In addition, image monitoring every 10 min using 
strong sky lines show no variation in TF tuning during the 
observations. 

For each night 100 bias and 100 well exposed dome flat- 
field images were taken (^30,000 ADU/pixel/image). Our 
program uses dome-flats, as opposed to sky-flats, as small 
scale pixel-to-pixel variations are a potentially important 
source of noise, requiring a large number of well exposed 
flats to effectively remove it. The dome lights are not very 
homogeneous, in wavelength or spatially, which complicates 
flat-fielding large-scale features, and an illumination correc- 
tion was added to the final co-added flat-field image. Our 
observations were all specified to occur at the same pixel lo- 
cation (no dithering), to help suppress flat-fielding errors, 
with the guiding each night keeping the target's location 
within a few pixels. 

2.3. Reduction 

The bias frames and flat fields were combined using stan- 
dard IRAF routines, and used to correct each image, be- 
fore aperture photometry was performed using IDL for the 
target and reference stars. The aperture location was de- 
termined by fitting a 2D Moffat function ( |Moffat|[l969l ) to 
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Fig. 2. GTC OSIRIS narrow-band transit light curves at (top-to-bottom) 8839, 7665, 7582, and 6792 A, with an arbitrary 
flux offset applied to each curve for clarity. The best-fit models are also over-plotted, along with the residuals (bottom four 
curves) which have a y-axis error bar corresponding to the per-point flux uncertainty. The light curves are near-Poisson 
limited, with per-point S/N levels between 1149 and 1351. 



the point spread function (PSF) of the target and refer- 
ence stars within each image. The background sky was 
subtracted using the pixels from an annulus around each 
aperture. After selecting a wide range of apertures and sky 
annulus ring sizes, aperture radius sizes between 18 and 20 
pixels were found to be optimal for our observations, along 
with sky annulus sizes in the ranges of 70 pixels for the 
inner sky ring radius and 150 pixels for outer radius. For 
a given night, the same sized apertures were used for both 
the target and reference star. From our observed PSFs, the 
theoretical maximum S /N aperture (which balances the in- 
creased readout noise of including more pixels with the de- 
creasing num ber of photons away from the PSF center, see 
Howell|2000 1 was found to be 20 pixels, in good agreement 
with the empirically determined optimal aperture size. For 
the Dec 2009 observations at 7582 A, the average count level 
was observed to be 3.44x 10 6 counts for the target XO-2 and 
3.55 x 10 6 counts for the reference star XO-2B. Including the 
errors from photon-noise, readout noise and sky subtrac- 
tion, the average per-exposure uncertainty is found to be 
cr XO2 =0.00053 for XO-2 and a REF =0.00053 for the refer- 
ence star, giving a differential photometric light-curve (LC) 
with an average per-exposure error of <tlc=0. 00074, or a 
S/N per image of 1351. The error including photon, read 
noise, and sky background uncertainties were tabulated for 
each exposure, and used as the estimate of the error in 
the differential photometric light-curve. Similar values are 
found for the other two nights, with an average S/N per im- 



Table 2. Three parameter non-linear limb darkening coef- 
ficients for GTC OSIRIS Tunable Filters 



Wavelength (A) 


6792 


7582 


7664.9 


8839 


C2 


1.4055 


1.4003 


1.4036 


1.4135 


C3 


-1.0570 


-1.3042 


-1.2957 


-1.3454 


C 4 


0.3453 


0.4802 


0.4745 


0.4958 



age of 1167 for the 6792 A exposures and an average S/N 
of 1342 for the 8839 A exposures. 

The potassium line-core has nearby strong telluric O2 
absorption lines, which could in principle affect the differ- 
ential photometry. However, we found that the resulting 
atmospheric differential extinction between the target and 
reference star was negligible, < 4 x 10~ 5 , when modeling 
the transmission function t for each sta r through the atmo- 
spher e as a function of airmass, £ sec ( z ) ( Vidal-Madjar et al. 
20T0b. 



3. Analysis 

3.1. Transit light curve fits 

We modeled the transit light curve w ith the analytical tran- 
sit models of |Mandel fc Agol| ( |2002[ ), fitting for the central 
transit time, planet-to-star radius contrast, limb-darkening, 
inclination, stellar density, and baseline flux. The errors on 
each datapoint were set to the tabulated values which in- 
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Table 3. System parameters for XO-2b from individual light curve fits 



Parameter 


Value 








central wavelength, A (A) 


6792.0 


7582.0 


7664.9 


8839.0 


wavelength range, AA(A) 


12 


12 


12 


12 


planet/star radius contrast, i?=R p ;/R* 


0.1043±0.0021 


0.1067±0.0015 


0.1083±0.0017 


0.1032±0.0012 


inclination, i [deg] 


88.4±1.8 


87.23±0.99 


86.09±0.79 


88.5±1.9 


impact parameter, 6=acosi/R* 


0.23±0.25 


0.37±0.13 


0.49±0.10 


0.21±0.21 


system scale, a/R, 


8.14±0.43 


7.69±0.36 


7.20±0.32 


8.05±0.35 


stellar density, p* [g cm -3 ] 


1.49±0.25 


1.26±0.17 


1.03±0.14 


1.44±0.24 


linear limb- darkening coeff., ci 


1.3943±0.048 


1.412±0.053 


1.491±0.061 


1.448±0.043 



Table 4. Joint fit system parameters for XO-2b 



Parameter 


Value 


Period, [days] 


2.61586178±0.00000075 


Mid-transit Time, T(0) f [HJD] 


2454466. 88457±0.00014 


radius reference, i? r e/ =R P i/R* 


0.10362±0.00070 


altitude, Z{\ = 6792) 


0.00146±0.00075 


Z(\ = 7582) 


0.00278±0.00077 


Z(X = 7664) 


0.00322±0.00079 


Z(\ = 8839) 


(fixed) 


inclination, i [deg] 


87.62±0.51 


system scale, a/R* 


7.83±0.17 


impact parameter, 6=acosi/R* 


0.324±0.070 


stellar density, p* [g cm -3 ] 


1.328±0.088 


linear limb- darkening coeff., ci 


fixed to Table |2| values 


T(261) [HJD] 


2455149. 62450T0. 00027 


T(269) [HJD] 


2455170.55140±0.00023 


T(269) [HJD] 


2455170.55181±0.00023 


T(303) [HJD] 


2455259.48978±0.00030 


Using best fit values (i, 


a/R*, b, & p*) 


R pi /R*(A)=.R(6792A) 


0.10508±0.00052 


i?(7582A) 


0.10640±0.00058 


i?(7664A) 


0.10686±0.00060 


i?(8839A) 


0.10362±0.00051 


Previously determined model dependent parameters 


R* [R©] 


0.976±0.020 :t 


Rpi [Rjup] 


0.996±0.025 t 


M* [M ] 


0.971±0.034 t 


M pl [M Jup ] 


0.565±0.054 t 



(2009) 



^fit along wi th the data from Burke et al. (|2007|) and Fernandez 
let al 



From Fernandez et al. (20091 



eluded photon noise, sky subtraction, and readout noise. 
To account for the effects of limb-darkening on the transit 
light curve, we adopted the three parameter limb-darkening 
law, 



m. 
i{±) 



= 1 - c 2 (l - n) - c 3 (l - M 3/2 ) - c 4 (l - V? 



(1) 



calculating the coefficients following Sing (2010) using the 
transmission through our narrow-band tunable filters (see 
table [2]). We used the three parameter non-linear limb- 
darkening law, as it best provides both the flexibility to cap- 
ture the inherent non-linearity of the limb darkening inten- 
sity profile and helps to avoid the d eficient ATLAS model 
values near the limb (see Sing||2010 for more details). This 
law has now successfully been used to model high preci- 
sion transit light curves in cluding data from HST/NICMOS 



flSinget al.|2009[ ), Spitzer ( |Hebrard et al 12010 1 
2010), and CoRot. In the fits, we also allowed the baseline 



Desert et al. 



flux level of each visit to vary in time linearly, described 
by two fit parameters. The linear trend accounts for pos- 
sible pixel-scale related drifts during the observations, and 
we found that two nights exhibited a small baseline gra- 
dient. The best fit parameters were determined simultane- 
ously with a Leven berg-Marquardt least-squares algorithm 
(Markwardt 2009) using the unbinned data. A few deviant 
points from each light curve were cut at the 3-cr level in 
the residuals, likely due to cosmic rays. The results from 
fitting each light curve separately are given in Table [3j In 
the separate light curve fits, we allowed the linear coeffi- 
cient term ci to vary freely, with the other two coefficients 
fixed to the predicted model values. Choosing to fit only one 
limb-darkening parameter ensures that the fits are both not 
significantly biased by adopting a stellar atmospheric model 
and do not suffer from degeneracies between fitting multi- 
ple limb-darkening parameters. The fit values for the linear 
limb-darkening coefficients are measured at the ~4% level, 
and are consistent within la of the model values (see Table 
[2]), indicating that the Atlas models and three parameters 
law are performing well. 

We also performed a joint fit for the four light curves, si- 
multaneously fitting for the inclination i and stellar density 
p*, with the limb-darkening values fixed to the predicted 
model values and the results plotted in Fig. [2j As the preci- 
sion and spectral resolution in our light-curves are sufficient 
to observe atmospheric features, which cause measurable 
changes in the planet's radii, we jointly fit each light curve 
with a common reference planetary radius, R re f=R p i/Tl+, 
as well an altitude in units of stellar radius, Z = Z(X), 
above the fit reference radius value with 



R{\) = Z(X) + R 



ref- 



(2) 



The transit light curve at 8839 A was used as the refer- 
ence radius value, setting the altitude Z(8839)=0, as it con- 
tained the smallest measured radius. The main advantage 
of jointly fitting R re f for the different color light curves, 
and separately fitting Z{\), is that the uncertainties in fit- 
ting for the system parameters (i, a/R*, 6, & p*) do not 
affect the measurement (joint probability distribution) of 
Z(X), given that the correlations between the system pa- 
rameters and radius are largely wavelength independent. 
We also quote the values for the individual R p i/R* values 
for each light curve, when the system parameters are set to 
their best fit values. The results of these joint fits are given 
in Table H 

3.2. Noise Estimation 

Each of our light curves achieve sub-mmag precision, with 
the RMS of the residuals between the best-fitting individual 
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models and the data being 0.087% at 6792 A (S/N of 1149), 
0.074% at 7582 A (S/N of 1351), 0.076% at 7664 A (S/N 
of 1315), and 0.083% at 8839 A (S/N of 1204). These val- 
ues compare very favorably to the theoretically estimated 
error in §2, dominated by photon noise, indicating that our 
transit light curves are regularly achieving between 90% to 
100% of the theoretical Poisson limit, without having to 
use de-correlation techniques to remove systematic errors. 
Other sources of noise are found to be negligible. This in- 
cludes scintillation, as was also found to be the case with 



the large aperture of GTC by |Colon et ah] < |2010 1. 

We checked for the pres ence of systematic errors corre- 
lated in time ("red noise", Pont et al. 2006) by checking 
that the binned residuals followed a iV -1 / 2 relation when 
binning in time by N points. In the presence of red noise, the 
variance can be well modeled to follow a a 2 = cr^/N + of 
relation, where a w is the uncorrelated white noise compo- 
nent while ay characterizes the red noise. We found no sig- 
nificant evidence for red noise, when binning on timescales 
up to the ingress and egress duration, with all the deter- 
mined ay values below 1 x 10 -4 , a level corresponding to the 
limiting precision that our radii values can be measured 
within an individual light curve (ct#!~1.1x 10~ 4 ). When 
possible, we also checked for the presence of red no ise with 
the wavelet-based method of Carter & Winn ( 2009 ), finding 
good general agreement with the binning technique. 



4. Discussion 

The system parameters determined with the GTC light 
curves are in good agreement (at the ^l-a level) with 
the most accurate values for XO-2b previously published 
in the literature, namel y the z-band measurements from 
Fernandez et al. (20091. We choose our wavelengths to 
be maximally sensitive to the potassium line, with the 
transits measured at 7664 and 7582 intended to measure 
the K I line core, while those at 6792 and 8834 A were 
chosen to have minimum contributive opacities from K, 
Na, and H2O. Our joint-fit system parameters measure a 
significant altitude z at two wavelengths, with the mea- 
sured potassium line center at 7664 A to be z(7664)= 
0.00322±0.00079 i?*, which is above the reference alti- 
tude at the 4.1-er significance level ("absorption" level of 
A(i? pi /i?^) 2 = 0.067±0.016%). Assuming the stellar ra- 



dius as given in Table [4j the potassium line core is seen to 
reach altitudes of 2190±540 km above the 8839 A reference 
planetary radius. A significant altitude is also measured at 
7582 A, with z(7582)=0.00278±0.00077 R* which is above 
the reference zero altitude at the 3.6-ct significance level. 

The wavelength-dependent transit radius is defined as 
the radius where the total slant optical depth reaches 
r eo =0 .56, as determined in Lecavelier des Etangs et al 
(2008) , with the transmission spectra models of Fortney 



et al. (2010) matching these results. For an isothermal at- 



mosphere, the effective altitude z of the atmosphere as a 
function of wavelength A is given in Lecavelier des Etangs 
et all (120081) to be 



z(A) 



m ( ea bs P ref <7(\) 2wRpi\ 
V kT ^9 J ' 



(3) 



' cq 




7000 

ilength (A) 



Fig. 3. XO-2b narrow-band GTC photometric results at 
four wavelengths ( colored points) a l ong w ith the z-band ra- 
dius (black dot) of |Fernandez et al.| ( [2009| . The X -axis error 
bars indicate the wavelength range of the photometry, while 
the Y-axis error bars indicates the errors on the measured 
planetary radius, assuming a stellar radius of 0.976 R . 
Also over-plotted is a 1500 K isothermal mo del (binned to 
R=50 0) of the transmission spectrum from Fortney et al. 
(2010) adapted for XO-2b, which is dominated by gaseous 
Na and K. TiO and VO opacities are removed in this model. 
The open circles indicate the predicted filter-integrated 
model values. 



is the absorption cross section, and H is the atmospheric 
scale height. The measured altitude parameter Z{\) from 
§3 relates to z(A) by Z(X)=z(X)/R ic . 

Assuming a temperature of 1500 K for XO-2, based 
on the broadband s econdary eclipse measurements of 
|Machalek et al. (2009), and a surface gravity of 14.71 m/s 2 , 
based on the planetary mass and radius values found by 
|Fernandez eTa l. (2009), we estimate an atmospheric scale 
height of #=370 km7and ff/.R*=0.00054. Thus, the in- 
dividual precisions of our wavelength dependent planetary 
radii measurements {(Jr i /r 1< ~0. 00055 & az ~ 0.00075) 
correspond to one atmospheric scale height. The atmo- 
spheric feature at 7664 A thus spans 6±1.5 scale heights 
above the reference altitude. 



where e a bs is the abundance of dominating absorbing 
species, P re f is the pressure at the reference altitude, u(A) 



4.1. Atmospheric signatures of potassium 

We compare our best-fit planetary radius values to the re- 
cent hot- Jup iter model atmospheric calculations of Fortney 
|et al.| p"010[ ), along with XO-2b specifically generated mod- 
els. These models include a self-consistent treatment of ra- 
diative transfer and chemical equilibrium of neutral species, 
though do not include the formation of hazes. Given our 
limited wavelength coverage, such physically motivated 
models are necessary for comparison with our data, as 
there is not enough data to justify freely fitting for a large 
number of atmospheric parameters. The comparisons here 
are an effort to more broadly distinguish XO-2b's atmo- 
sphere between different classes of models, given there are 
no prior constraints on the basic composition or presence of 
clouds/hazes, as well as to compare the specific predictions 
these models have for XO-2b. 
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TiO/VO 




7000 
Wavelength (A) 



Fig. 4. Same as Fig. |3j but over-plotted with a 1500 K 
isothermal model also containing TiO/VO opacity, produc- 
ing many strong absorption features. The TiO/VO abun- 
dances in the model are from equilibrium chemistry. There 
is no evidence for TiO/VO in the GTC data, as the model 
fits with TiO/VO are significantly worse than the alkali- 
dominated models (see text). 



The transmission spectrum calculations are performed 
for ID atmospheric pressure-temperature (P-T) profiles. 
The model atmospheres use the equilibrium chemical 



abundances, at sol ar metallicity, described in |Lodders ~fc 
Fegley fl2002 2006[). Th e opacity database is described in 
Freedman et ai. ([20081. The model atmosphere code has 
been us ed to generate P-T profiles for a variety of close-in 
planets ( |Fortney et al.|2005|[2006l|2008| . The transmission 
' ' ' Fc 



spectrum calculations are described in |Fortney et al. (2010). 
For XO-2b, we find that the optical opacity is dominated by 
gaseous Na and K, as is expected for many close-in planets 
over a broad temperatu re range (e.g., |Sudarsky et al.|2000 



Seager & Sassclov 2000 ) . These alkalis may be subjected to 
photo-ionization, but this is not included in this study. 

The models are calculated wavelength-by- wavelength at 
a resolution of R=3500, about seven times higher than the 
resolution of our data. For the comparisons to the GTC 
data, we assumed a value of i?*=0.976 M Q (Table 3), and 
calculated the filter-integrated predictions from the high 
resolution model, fitting only for an arbitrary model alti- 
tude shift. Our GTC measurements are well matched with 
theoretical models which include significant opacity from 
atmospheric potassium. The first indications to date that 
this species is present in hot- Jupiter atmospheres. 

Isothermal models: Comparing the observed atmo- 
spheric features to isothermal hydrostatic uniform abun- 
dance models helps provide an overall understanding of the 
features observed, and departures from those conditions. 
Using equation 3, we adapted generic hot-Jupiter isother- 
mal models to values appropriate for XO-2. We find an 
isothermal model dominated by alkali metals in the optical 
is a substantially better fit than a grey "flat" atmosphere, 
as would be expected if high-altitude clouds were dominant. 
Adjusting only the model reference altitude (i.e. the model 
normalization) in the fits to the GTC and z-band points 
gives a \ 2 01 22 with 4 degrees of freedom (DOF) for a 
grey featureless atmosphere and a x 2 of 7.0 for a 1500 K 




700 800 900 1000 1100 1200 1300 1400 1500 1600 1700 1800 1900 
Temperature (K) 

Fig. 5. Temperature-Pressure profiles for the XO-2b spe- 
cific models, including a planet-wide averaged profile (left) 
and a dayside averaged profile (right). The GTC transit 
data are sensitive to pressures from ~10 _1 bar at the low- 
est altitudes to ^10~ 3 bar at the highest altitudes in the 
potassium and sodium line-cores. 



isothermal alkali-dominated atmosphere (see Fig. [31 with 
the model binned to R=500). With a smaller scale neight, 
a cooler 1000 K isothermal model somewhat under-predicts 
the K I core levels, producing a % 2 fit of 9.3 with 4 DOF. 

Models including TiO and VO in addition to K I do not 
improve the fits. A 1500 K isothermal model including TiO 
and VO opacity at chemical equilibrium abundances gives 
a x 2 n t of 12.6 with 4 DOF (see Fig. WJ). In these models, 
TiO /VO effectively raise the altitudes of the non-potassium 
core model levels, degrading the fit. These models still pre- 
dict the 7664 A feature due to potassium, as the low levels 
of TiO do not cover the potassium core signature, with the 
K-doublet still visible in Fig. [4] Models completely domi- 
nated by only TiO/VO absorption, perhaps appropriate to 
very hot-Jupiters such as Wasp-12, can be definitively ruled 
out, as they provide a substantially poorer fit to the data 
(x 2 of 32). Sulfur compounds are unlikely to contribute at 
these longer wavelengths. 

XO-2b specific models: We also compared our GTC 
planetary radius values to models specifically calculated for 
the known parameters of XO-2b, assuming chemical equi- 
librium at solar metallicity. A model using a dayside av- 
eraged temperature-pressure (T-P) profile (see Fig. [5] and 
Fig. K3J) is our best fitting model, giving a % 2 value of 6.6 
for 4DOF. The cooler model transmission spectra using 
a planet- wide averaged T-P profile (see Fig. [5]), performs 
moderately worse giving a fit with a x 2 of 8.0 for 4 DOF. 
The planet- wide averaged model predicts lower overall tem- 
peratures than the day-side model, also resulting in a slight 
under-prediction of the potassium line-core signature simi- 
lar the 1000 K isothermal model. 

Both XO-2b specific models and the non-TiO isothermal 
models predict that the opacity from potassium, sodium, 
and H2 Rayleigh scattering dominate the optical, with H2O 
absorption features prominent red-ward of ^9000 A. From 
our data, the observed strong potassium opacity in the op- 
tical and lack of significant TiO is in general agreement 
with XO-2b being classified as a transition "pM/pL" hot- 
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Fig. 6. Same as Fig. [3j but over-plotted with a XO-2b specific model of the transmission spectrum, which uses a 
dayside-average temperature-pressure profile. The model is dominated in the optical by Na and K absorption signa- 
tures. Compared to the isothermal model of Fig. [3J the lower temperatures at higher-altitudes in this model lead to 
a better fit of the K-core 7664.9 A measurement. Further observations are planned at blue wavelengths (blue boxes), 
which will be highly sensitive to atmospheric sodium, and will help provide a more complete picture of the transmission 
spectrum. 



Jupiter. While the GTC data does not provide particularly 
strong temperature constraints, we find that temperatures 
in the range of ^1500 K from the isothermal and day-side 
models offer slightly better fits than the cooler planet-wide 
and 1000 K models. Using 1500 K temperatures or a specific 
T-P profile from the dayside would at first seem at odds 
with transit spectral data, associated with the planetary 
limb. However, 3D hot- Jupiter models have recently shown 
that the warmer dayside temperatures can increase the at- 
mospheric scale height and effectively "puff up" the dayside 
atmosphere, covering both the cooler plane tary limb and 



night -side transmission spectral signatures (Fortney et al. 



20101. Our optical transit data are sensitive to pressures 
from ^10 _1 bar at the lowest altitudes between the Na I 
and K I line wings down to ~10 -3 bar at the highest al- 
titudes in the potassium and sodium line-cores (assuming 
solar abundances) . This is roughly the same pressure range 
that Spitzer secondary-eclipse emission spectra measure- 
ments are sensitive to, given that the smaller cross sec- 
tion in the optical (compared to the IR) is compensated by 
the larger horizontal integrated density of the slant-transit 
geometry. These similarities have the potential to help in 
identifying unknown absorbers thought to be responsible 
for day-side inversions, as thei r signatures cou l d app ear 



ported a "weak" temperature inversion layer for XO-2b, 
with models placing it near a 10 -2 bar level. Our models 
lack thermal inversions, though we do find similar ^1500 
K temperatures as models containing a thermal inversion 



( Machalek et al.|2009 ). Given that non- inverted models still 
fit the Spitzer IRAC photometry at 2.2-ct levels or better, 
and the average limb profile may still be sufficiently differ- 
ent than the day-side and lack detectable inversions, we do 
not consider the difference as significant. Further secondary 
eclipse measurements will be necessary to definitively de- 
termine the presence or absence of a thermal inversion layer 
in XO-2b. 



4.2. Limits on stellar activity 

A transmission spectrum in active stars can be affected 
from the presence of stellar spots, both when a spot is oc- 
culted by the planet and by the presence of non-occulted 



et al. 


2008 


Agol et al. 


2010) 



in optical transmission spectra. Machalek et al. (2009) re 



is normalized, such that different epochs with different 
amounts of spots will show small but non-negligible radius 
changes. For the XO-2b observations, sufficiently large stel- 
lar activity could in principle affect the measured radii at 
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the reference wavelengths (6792 & 8839 A) compared to 
the potassium line core wavelengths, as they were not taken 
simultaneously during a transit event. However, in the ab- 
sence of strong stellar activity, multiple transits spanning 
long epochs (i.e. non-simultaneous transits) can still be used 
and fit together at high precision, without the need for a 
stellar activity related correction. 



Czesla et al. (120091) and Desert et al. (2010) proposed 
that the impact of rotation modulated spots on the final 
fitted radius could be estimated following 



ARpi/R± 

r t 



<*f\ 

2 : 



(4) 



where fx is the relative variation of the stellar flux due 
to stellar activity (i.e. rotation modulated spots) and a 
is an empirically determined constant of order unity. In 
the hypothesis that the stellar surface brightness outside 
dark stellar spots are not modified by the activity level, 
the parameter a is given by —1. The range of fx between 
active and non-active stars is very large, spanning many 
orders of magnitude. For the highly active star HP 189733, 
fx ~ 2 - 3% in the visible ( |Henry fe W^[2008| ), while 
for the non-active star TrES-2, fx was found to have an 
amazingly low value of (5 . 18±0.19) x 10~ 6 from Kepler data 
( |Kipping fe Bakos||2010| . |Ciardi et ah] ( |2010[ ) studied the 
stellar photometric variability for stars in the Kepler sam- 
ple, finding a bi-model distribution between active an non- 
active stars. K-dwarfs with low activity (2/3 of the sample) 
were found to typically display a stochastic photometric 
variability of less than 1 mmag, while active stars were of- 
ten found to be periodic on timescales of days to weeks. 
For the HS T transit observations of the high l y active star 



HP 189733 (Pont et al. 2008 Sing et al. 2009 Desert et al. 



20101, an observational strategy was adopted to counter 



this effect, which involved monitoring the stellar activity 
level. In this way, the stellar baseline flux was known at 
every epoch of the transit observations, and could then be 
used to correctly normalize the data. For non-active stars 
such as TrES-2 or XO-2, this strategy may not be needed 
as the transit-fit planetary radii values will not significantly 
change from one epoch to another. 

With the TrES-2b Kepler data, and such a low level 
of activity, eighteen opti cal transits were fit together by 
Kipping & Bakos (2010), who found no signs of epoch- 
dependent radius changes, as would be apparent from 
strong activity-related effects. The star TrES-2 has very 
similar activity indices as XO-2, with |Knutson et aT] ( |2010 ) 
measuring the stellar activity indices of the Ca 11 H & K 
line strength to be log(i?^)=-4.988 for XO-2 and -4.949 
for TrES-2, while the S- value of XO-2 is 0.173 and 0.165 for 
TrES-2. Similar activity levels as TrES-2 suggests activity- 
related radii changes in XO-2 are also negligible. The fi- 
nal TrES-2 combined transit light curve was successfully fit 
down to the 20 ppm level, which is nearly an order of mag- 
nitude below our required levels for XO-2b (~100 ppm). 
Hat-P7 and Kepler-4,5,6,7,8 also show a low photometric 
variability around the ^0.2 mmag level, with the t ransits 
removed wh en calculating the light curve statistics ( Ciardi 
et al. 2010[), an d all have similar activity indices as XO-2 



QKnutson et~aT|2010[ ) 



For XO-2, there is also no obvious evidence for occulted 
spots in any of our sub-mmag precision light curves con- 
sistent with a non-active star, though non-occulted spot 



features could be clustered at different latitudes. In addi- 
tion, no spot-induced rotational modulations have yet been 
reported, with the discovery XO-2 light-curves see mingly 
constant do wn to a noise- limited ^7 mmag level (Burke 
eT^|2007] . Assuming fx < 0.007 for XO-2 and a = -1, 
we find A(R p i/R+) < 0.00036, which is below the limiting- 
precision level of our measurements, and below the esti- 
mated size of one atmospheric scale height. While current 
evidence points towards XO-2 as being a non-active star 
with a likely small photometric variability, further activity 
monitoring of this and other transit hosting stars can play 
an important role in assessing the true variability. 

4.3. Transit Timing Studies 

We used the four central transit times of our GTC data 



along with the transi t times of Burke et al. (2007) and 
Fernandez et al. (2009), to determine the transit ephemeris, 
using a linear function of the period P and transit epoch 
E, 



T(E) = T(0) + EP. 



(5) 



The results are given in Table [4j where we find a period 
of P = 2.61586178 ± 0.00000075 (days) and a mid-transit 
time of T{0) = 2454466.88457±0.00014 (HJP). A fit with 
a linear ephemeris is in general an unsatisfactory fit to the 
data, having a \ 2 of 42 for 20 deg rees of freedom (\f, =2.1). 
The fit period is within l-er of the Burke et al. ( 2007 ) value 
of 2.615857±0.000005 days but 3-a away fronilne period 
of 2.615819±0.000014 days |Fernandez etal] ( [2009] ) found 
when fitting for the ephemeris wi th just their six transi t 
times. Similar to the conclusions of Fernandez et al. (2009), 
we find it is not clear whether to interpret the poor linear fit 
and differences in period to a genuine period variation from 
the presence of a second planet, or whether the uncertain- 
ties in all the timings are underestimated by ~ \/2. Given 
that the observed variations are so close to the per-transit 
noise level, an unambiguous transit timing variation and 
detection of a second planet free of systematic errors will 
require more observations with generally improved timing 
errors. GTC OSIRIS observations in sub-array mode should 
help provide the necessary precision, as the duty cycle can 
be increased to near 100%. 



5. Conclusions 

We present GTC OSIRIS observations for three transits of 
XO-2b, at four carefully selected narrowband wavelengths, 
which have revealed significant absorption by atmospheric 
potassium. This is the first evidence for this important al- 
kali metal in an extrasolar planet, which has long been theo- 
rized along with sodium to be a dominant source of opacity 
at optical wavelengths. Identification of atmospheric con- 
tents is one of the first steps to understanding the nature 
of exoplanetary atmospheres. The presence (or lack) of key 
species can give clues on the exoplanet's temperature, pres- 
ence of condensate clouds, chemistry, and can potentially 
be used to type different atmospheres and place exoplan- 
ets within different sub-classes. Clear gaseous hot- Jupiter 
atmospheres are thought to have low geometric albedos, 
appearing very dark at visible wavelengths, with broad 
pressure broadened potassium and sodium line- wings effi- 
ciently absorbing the incoming stellar radiation (Sudarsky 
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et al.|2000| Seager et al.|2000 ). A detection of potassium in 
XO-2b helps support this overall picture. Potassium and 
sodium could also play important roles in infla ting hot 
Jupiters through ohmic dissipation mechanisms (Batygin 
& StevensonpOlO). Notably, there is also no particular ev- 



idence for K I photoionization, which could have been evi- 
dent by depletion at high altitudes. 

Narrow-band spectrophotometry is a robust method for 
obtaining precision differential photometric transit light- 
curves which have a spectral resolution sufficient to detect 
atmospheric transmission spectral features. OSIRIS is the 
first instrument capable of producing sub-mmag high-speed 
spectrophotometry, and coupled with the large 10.4 meter 
aperture of GTC, is capable of detecting spectral features in 
transmission spectra on exoplanets such as XO-2b, whose 
host star is 3.5 magnitudes fainter in the optical than those 
of the well studied transiting hot Jupiters HD 189733b and 
HD 209458b. XO-2b now joins those two exoplanets as 
the only hot- Jupiters where direct constraints on the at- 
mospheric opacity in the optical can be made. Detecting 
and studying the atmospheres for such fainter targets will 
be critical for a broader understanding of the hot- Jupiter 
class as a whole, as most of the known transiting planets 
have been found at magnitudes >10. Further observations 
from our program, including planned XO-2b measurements 
at blue wavelengths, can help conclude whether the opti- 
cal atmospheres of these exoplanets are dominated by the 
two alkali metals sodium and potassium together, search 
for Rayleigh scattering, and place limits on other optical 
absorbers thought to be linked to thermal inversion layers. 
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Note Added. Shortly before submission, we became 
aware of i ndependent, but si milar observations of another 
exoplanet ( Colon et al.||2010 ). Future observations of these 
and additional exoplanets will enable comparisons of the 
atmospheric composition and structure, as well as stud- 
ies of potential correlations with other planet or host star 
properties. 
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